ON  THE  ACCURACY  ANALYSIS 
OF  AIRBORNE  TECHNIQUES 
FOR  PASSIVELY  LOCATING 
ELECTROMAGNETIC  EMITTERS 

L.  H.  Wegner 


A  Report  prepared  for 

UNITED  STATES  AIR  FORCE  PROJECT  RAND 

I<preduwl  toy 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

IprinsflaVd.  V».  MIS1 


Rand 

SANTA  MONICA  CA  90406 


7? 


DOCUMENT  CONTROL  DATA 


I.  ORIGINATING  ACTIVITY 


The  Rand  Corporation 


2a.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2k.  GROUP 


ON  THE  ACCURACY  ANALYSIS  OF  AIRBORNE  TECHNIQUES  FOR  PASSIVELY  LOCATING 
ELECTROMAGNETIC  EMITTERS 


4.  AUTHORS)  (last  mm,  f Inf  nama,  Initial) 

Wegner,  L.  H. 


.  ABSTRACT 


Discussion  o£  methods  of  locating  electro¬ 
magnetic  emitters  from  airborne  electronic 
reconnaissance  systems,  using  passive  mea¬ 
surements  on  the  electromagnetic  waves  of 
the  emitter  taken  from  several  locations. 
Measurements  for  finding  the  emitter  may 
be  either  the  direction  of  arrival  at  2 
or  more  locations  along  the  aircraft's 
flight  path  for  single- aircraft  DF  (direc¬ 
tion-finding)  systems,  or  the  relative 
time  of  arrival  of  emitter  pulses  at  the 
different  aircraft  locations  for  3-aircraft 
TOA  (time-of-arrival)  c vs terns.  Emitter 
location  accuracy  is  a  function  of  the 
accuracy  of  these  indirect  measurements 
of  the  emitter  location  and  the  estimation 
procedure  used  to  combine  the  accuracies 
of  various  measurements  of  location  and 
bearing.  A  general  expression  is  derived 
for  a  lower  bound  to  the  covariance  matrix 
of  the  unbiased  estimates  of  emitter  lo¬ 
cation  coordinates  (and  consequently  the 
location  CEP),  using  the  Cramer-Rao  ine¬ 
quality.  Examples  illustrate  the  applica¬ 
tion  of  the  Cramer-Rao  lower  bound  to  the 
emitter  CEP  for  both  DF  and  TOA  systems. 
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PREFACE 


A  primary  function  of  electronic  reconnaissance  is  to  detect  and 
locate  enemy  radars  and  other  electromagnetic  emitters  such  as  those  in 
communications  networks.  Knowledge  of  the  location  of  radars  in  enemy 
air  defense  systems  is  necessary  for  attacking  the  Individual  radar 
sites  and  associated  equipment  such  as  surface-to-air  missiles  and  for 
developing  effective  aircraft  penetration  tactics.  In  the  battlefield 
area  the  location  of  enemy  radars  and  communications  emitters  can  be 
used  to  locate  individual  ground  units  ac  well  as  to  provide  important 
inputs  to  the  overall  intelligence  picture. 

As  an  aid  in  comparing  the  relative  effectiveness  of  alternative 
airborne  emitter  location  finding  systems,  a  computer  model  was  devel¬ 
oped  at  The  Rand  Corporation  to  simulate  the  capabilities  of  such  sys¬ 
tems  against  postulated  full-scale  emitter  environments.  An  important 
output  of  the  model  is  the  accuracy  with  which  each  detected  emitter  is 
located.  This  report  provides  the  documentation  for  the  methods  used 
in  the  model  for  determining  emitter  location  accuracy. 

This  report  should  be  of  interest  to  individuals  involved  in  ana¬ 
lyzing  and  evaluating  emitter  location  finding  systems.  It  presents 
general  methods  for  determining  location  accuracy  and  provides  a  number 
of  analytic  and  numerical  results  which  facilitate  an  understanding  of 
the  effects  of  system  component  measurement  errors  and  aircraft/emitter 
geometry  on  overall  system  location  accuracy.  Although  most  of  the 
mathematical  derivation*:  are  Included  for  completeness  and  to  present 
a  unified  treatment  of  the  location  accuracy  problem,  the  reader  who  is 
not  interested  in  the  mathematical  development  can  make  practical  use 
of  the  results  and  curves  presented  to  make  his  own  rapid  determina¬ 
tions  of  emitter  location  CEP. 


SUMMARY 


This  report  treats  methods  of  locating  electromagnetic  emitters 
from  airborne  electronic  reconnaissance  systems,  using  passive  mea¬ 
surements  on  the  electromagnetic  vaves  of  the  emitter,  taken  from 
several  different  locatione.  The  measurements  used  for  location  find¬ 
ing  may  be  either  the  direction  of  arrival  at  two  or  more  locations 
along  the  aircraft's  flight  path  for  single-aircraft  DF  (direction- 
finding)  systems,  or  the  relative  time  of  arrival  of  emitter  pulses 
at  the  different  aircraft  locations  for  three-aircraft  TOA  (time-of- 
arrival)  systems. 

Emitter  location  accuracy  is  a  function  of  the  accuracy  of  these 
"indirect"  measurements  on  the  emitter  location,  the  accuracy  of  the 
aircraft  location  measurements,  the  accuracy  of  tha  bearing  measurement 
base  line  (in  the  case  of  DF  systems),  and  the  estimation  procedure 
used  to  combine  the  various  measurements.  A  general  expression  for  a 
lower  bound  to  the  covariance  matrix  of  the  unbiased  estimates  of  the 
emitter  location  coordinates  (and  consequently  the  location  CEF)  is  de¬ 
rived,  using  the  Cramdr-Rao  inequal icy.  The  lower  bound  is  independent 
of  the  estimation  procedure  used;  however,  it  is  shown  (in  an  appendix) 
that  the  covariance  matrix  of  the  "generalized  least  squares"  estimate 
of  the  location  coordinates  approximates  the  lower  bound.  Thus,  the 
lower  bound  may  be  used  as  an  intrinsic  measure  of  the  location  accu¬ 
racy  possible  under  the  given  measurement  conditions. 

Examples  are  presented  illustrating  the  application  of  the  Cramdr- 
Rao  lower  bound  to  the  emitter  location  estimate  CEP  for  both  DF  and 
TOA  systems.  Approximate  analytic  equations  for  the  emitter  CEP  are 
derived  for  both  DF  and  TOA  systems  for  the  case  where  the  aircraft 
location  errors  are  negligible  so  that  the  emitter  location  errors 
are  due  solely  to  errors  in  the  measurements  on  the  emitter.  For  DF 
systems,  analytic  results  are  also  derived  for  the  case  where  measure¬ 
ments  on  the  elevation  angle  from  the  aircraft  to  the  emitter  are  taken 
in  addition  to  bearing  measurements.  For  TOA  systems,  additional  ana¬ 
lytic  results  are  presented  for  two  different  aircraft  location  schemes 
in  which  aircraft  location  accuracy  combines  in  a  simple  manner  with 
time-of- arrival  accuracy  in  determining  overall  emitter  location  accuracy. 
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lli  each  instance  where  analytic  results  are  obtained,  data  are  pre¬ 
sented  in  the  form  of  graphs  and  simplified  equations  which  enable  the 
user  who  is  not  concerned  with  the  background  mathematics  to  quickly 
estimate  emitter  location  CEPs  for  a  number  of  cases  of  practical  in¬ 
terest.  For  both  DF  systems  and  TOA  systems,  the  general  procedure 
using  the  Cramdr-Rao  lower  bound  to  determine  an  approximate  emitter 
location  CEP  from  the  combined  accuracies  in  the  emitter  measurements 
and  the  aircraft  location  measurements  is  illustrated  by  sample  cal¬ 
culations. 
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I.  INTRODUCTION 


The  problem  of  locating  an  object  on  the  surface  of  the  earth  by 
the  use  of  Indirect  measurements  has  a  long  history.  Surveyors  and 
navigators*  for  example,  have  always  been  concerned  with  location  as** 
timatlon.  Analytical  methods  for  the  solution  of  such  problems  can 
be  tracad  back  to  Causa  and  Laplace. ^  This  report  is  concerned  with 
the  analysis  of  the  accuracy  with  which  surface-baaed  electromagnetic 
emitters  can  be  located  by  passive  measurements  on  electromagnetic 
waves  from  the  emitter  taken  from  single-  and  multiple-aircraft  elec¬ 
tronic  reconnaissance  systems. 

To  locate  an  object  on  the  surface  of  the  earth  in  a  given  three- 
dimensional  coordinate  system  requires  measurements  on  at  least  three 
independent  quantities  functionally  related  to  the  throe  coordinate 
values  of  the  object.  For  locating  electromagnetic  emitters  passively* 
the  measurements  generally  taken  ate  emitter  altitude  (from  local  topo¬ 
graphic  data)  and  either  of  the  following: 

o  Direction  of  arrival  of  the  electromagnetic  waves  at  two  or 
more  locations  of  a  single  aircraft 

o  Time  of  arrival  of  emitter  pulses  at  three  aircraft  locations 

These  measurements  are  than  coablned  with  the  known  (or  estimated)  lo¬ 
cations  from  which  they  were  taken*  the  functional  relations  among  the 
various  locations  and  measured  quantities*  and  assumptions  about  the 
error  distributions  of  the  measurements,  to  arrive  at  an  emitter  loca¬ 
tion  estimate. 

When  the  emitter  measurements  are  taken  from  aircraft,  the  coor¬ 
dinates  of  the  aircraft  locations  as  well  as  the  emitter  location  must 
either  be  known  or  else  they  must  also  be  estimated.  For  each  unknown 
aircraft  location,  measurements  on  a  minimum  of  three  additional  quan¬ 
tities  functionally  related  to  the  aircraft  location  are  required. 

The  Air  Force's  ASQ-96  and  QRC-334  systems  are  examples  of  two 
basically  different  emitter  location  finding  systems.  The  ASQ-96  is 
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of -arrival  (directional  bearing)  measurements  on  emitter  pulses  from 
two  or  more  locations  along  the  aircraft’s  flight  path.  (If  desired, 
the  bearing  measurement  data  may  be  combined  with  those  of  a  second  DF 
aircraft  via  a  data  link  in  order  to  obtain  near-instantaneous  location 
estimates.)  The  aircraft  location  and  the  bearing  measurement  base 
line  are  provided  by  a  Loran-inertial  navigation  system.  In  the  QRC-334 
system,  time-of-arrival  (TOA)  measurements  are  made  on  emitter  pulses 
from  three  aircraft.  The  aircraft  locations  are  determined  from  mea¬ 
surements  of  the  ranges  between  aircraft  and  to  two  ground  stations  by 
DME  (distance  measuring  equipment)  systems,  and  aircraft  altitudes  from 
altimeters . 

In  the  "just-determined"  case,  in  which  there  are  exactly  as  many 
functionally  independent  quantities  measured  as  there  are  coordinates 
to  estimate,  the  location  estimation  problem  is  simply  one  of  determin¬ 
ing  the  (generally)  unique  solution  of  the  functional  equations.  How¬ 
ever,  when  there  are  more  measurements  than  quantities  to  estimate,  the 
"overdetermined"  case,  the  question  arises  as  to  the  proper  method  of 
combining  the  various  measurements  to  obtain  the  most  accurate  location 
estimate.  Solutions  for  this  problem  depend  upon  the  choice  of  crite¬ 
rion  for  the  average  "closeness"  of  an  estimate  as  well  as  the  probabil¬ 
ity  distributions  of  the  measurement  errors  and  the  relations  among  the 
measured  quantities. 

One  standard  measure  of  the  "closeness"  of  a  location  estimation 
procedure  is  the  location  CEP  (circular  error  probable),  the  circle 
around  the  true  location  within  which  50  percent  of  such  estimates  would 
lie.  When  the  emitter  coordinate  estimates  have  a  multivariate  normal 
distribution  with  mean  values  equal  to  the  true  coordinates,  the  loca¬ 
tion  CEP  can  be  expressed  as  a  function  of  the  variances  and  covariances 
of  the  location  coordinate  estimates.  In  Section  II,  a  lower  bound  is 
given  for  the  location  estimate  covariance  matrix  for  unbiased  location 
estimates.  The  lower  bound  is  derived  by  use  of  the  Cram4r-Rao  inequal¬ 
ity  under  the  assumption  that  the  measurements  have  a  multivariate  nor¬ 
mal  distribution  with  mean  values  equal  to  the  quantities  measured. 

In  Appendix  A  it  is  shown  that  the  covariance  matrix  for  generalized 
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least  squares  location  estlnates  approximates  that  given  by  the  lower 
bound.  In  Sections  III  and  IV,  emitter  location  CEPs  are  obtained  for 
DF  systems  and  TOA  systems  using  the  location  estimate  covariance  ma¬ 
trices  obtained  from  the  Cram£r-Rao  lower  bounds  for  the  system  con¬ 
figuration  examined. 
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II.  THE  ANALYTIC  BACKGROUND 

In  Che  most  general  form  considered  in  this  report,  the  emitter 
location  estimation  problem  is  as  follows.  Measurements  are  made  on 
an  unknown  emitter  location  (xq,  yQ,  zq)  from  N  aircraft  locations 
^xl’  yi‘  zi^'  *  *  * »  (*].  yN.  *N).  A  total  of  n  measurements.  Including 
the  emitter  measurements,  are  made  on  variables  functionally  related 
to  the  p  ■  3(N  +  1)  location  coordinates.  The  measurements  may  be 
either  "direct,"  i.e.,  measurements  on  the  individual  location  coor¬ 
dinates,  or  "indirect,"  i.e.,  related  functionally  to  the  location 
coordinates.  The  p  location  coordinates  may  be  considered  to  be  func¬ 
tionally  independent  parameters  to  be  estimated  and  the  n  measured 
variables  to  be  dependent  parameters.  In  general,  to  estimate  the  p 
location  coordinates  requires  at  least  p  relevant  measurements. 

For  a  convenient  notation,  let  the  vector 

0  -  c^,  *••,  ep)* 

-  <v  v  zo'  v  v  V* 

be  the  parameter  vector  to  be  estimated,  i.e.,  the  location  coordinates; 
the  vector 

m  -  (n^,  mn)* 

be  the  measurement  vector  on  the  parameter  vector 

P  “  (vlt  •••,  Pn)* 

i.e.,  m  is  the  measured  value  of  p  with  error  m  -  p;  and 

^(p.  8)  ■  0  i  -  1,  ••*,  r 

•f* 

All  symbolic  vectors  in  this  report  are  column  vectors,  so  the 
transpose  of  a  vector  v,  denoted  by  v  ,  is  a  row  vector. 


with  vector  form 


F(m »  6)  -  (fjCw,  B),  •••,  fr(p,  6))*  -  0  (1) 

be  the  functional  relationship*  among  the  components  of  p  and  0. 

Although  the  ultimate  goal  1*  to  obtain  estimates  of  the  emitter 
location  (xq,  yQ,  *0),+  it  it  convenient  to  consider  the  more  general 
problem  of  estimating  all  of  the  components  of  B.  If  there  are  direct 
measurements  on  the  aircraft  coordinates,  i.e.  ,  -  0^  for  some  i  and 

j,  then  the  corresponding  0j  may  be  omitted  from  B,  if  desired.  When 
the  estimates  of  the  aircraft  locations  are  also  of  interest  and  there 
are  direct  measurements  on  the  aircraft  locations,  it  may  be  advanta¬ 
geous  to  consider  the  more  general  problem  of  estimating  both  p  and  6 . 
This  is  done  in  the  Appendices. 

LEAST  SQUARES  AMD  MAXIMUM  LIKELIHOOD  ESTIMATES^ 

One  general  estimation  procedure  is  that  of  "ordinary  least 
squares in  which  the  estimates  p  and  B  are  the  values  of  p  and  0 
which  minimize  the  quadratic  form 

* 

■  (m  -  p)  (m  -  u) 

subject  to  the  constraints  F(y,  6)  ■  0. 

When  the  measurements  are  random  variables  with  mean  vector  p  and 
a  known  covariance  matrix 

I  -  E(m  -  p)  (m  -  p)* 


A  A 

the  generalized  least  squares"  estimates  p,  0  are  the  values  of  p  and 
B  which  minimize  the  quadratic  form 


ft  —1 

Q2  -  (m  -  p)  l  (m  -  p) 


(2) 


In  the  remainder  of  this  report,  ve  shall  generally  take  z  to  be 
the  altitude  coordinate  and  assume  sQ  to  be  known  with  negligible  error. 
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subject  to  the  constraints  F(p,  0)  -  0.  Generalized  least  squares  es- 
timafps  will  usually  be  better  estimates  than  ordinary  least  squares 
estimates. 

It  will  be  assumed  throughout  this  report  that  the  measurement 
vector  m  is  a  random  vector  with  a  multivariate  normal  distribution 
with  mean  vector  p  and  known*  covariance  matrix  E.  The  probability 
density  function  of  m  is  then 

♦  (»)  -  |Ef^(2n)"n/2  exp  [-$<»-  p)Vl(m  -  p>]  (3) 

"Maximum  likelihood"  estimates  of  p  and  S  are  the  values  of  p  and 
B  which  maximize  $(m)  subject  to  the  constraints  F(p,  B)  ■  0.  From 
the  form  of  <Km) ,  it  is  easy  to  see  that  maximizing  $(m)  is  equivalent 
to  minimizing  the  quadratic  form  in  the  exponent  of  (3) , 

Q3  *  (m  -  p)*2_1(m  -  p) 

subject  to  the  given  constraints,  F(p,  0)  -  0.  Since  Q3  is  identical 
to  Q2,  maximum  likelihood  estimates  are  identical  to  generalized  least 
squares  estimates  when  the  measurement  vector  has  a  multivariate  nor¬ 
mal  distribution.  This  will  not  be  true  in  general  for  other  proba¬ 
bility  distributions.  Maximum  likelihood  estimates  have  the  desirable 
property  that,  under  rather  general  conditions,  they  yield  estimates 
which  have  the  smallest  variance  when  large  numbers  of  measurements 
are  taken. 

The  determination  of  generalized  leaBt  squares  estimates  of  p  and 
B  requires  finding  the  minimum  of  the  quadratic  form,  Eq.  (2),  subject 
to  (possibly)  nonlinear  constraints,  F(p,  B)  *  0.  In  general,  this  re¬ 
quires  numerical  methods — one  general  procedure  is  described  in  Ap¬ 
pendix  A.  Since  the  primary  concern  of  this  report  is  with  emitter 
location  accuracy  and  not  estimation  procedures  per  se,  most  of  the 

1* 

In  practice,  I  is  obtained  from  error  analyses  and  independent 
testing  of  the  location  finding  system  and  its  components;  or  from 
hypothetical  values  when  different  proposed  systems  are  being  compared 
in  analytic  studies. 


remainder  of  the  discussion  Is  concerned  with  the  determination  of  lo¬ 
cation  accuracy. 


THE  CRAMER- RAO  INEQUALITY 

Maximum  likelihood  and  ordinary  and  general  least  squares  are  par¬ 
ticular  procedures  for  obtaining  estimates.  Other  estimation  procedures 
exist  and  r.d  hoc  estimates  not  based  on  any  general  procedure  are  used 
in  many  problems.  Various  criteria  are  available  for  judging  the  "good¬ 
ness"  of  an  estimation  procedure.  These  usually  involve  some  measure 
of  average  closeness  of  the  estimate  to  the  quantity  estimated.  For 
unbiased  estimates,  one  commonly  used  measure  of  average  closeness  is 
the  variance  of  the  estimate.  Under  fairly  general  assumptions  the 
Cramdr-Rao  inequality  described  below  provides  a  useful  lower  bound  to 
the  covariance  matrices  of  a  large  class  of  unbiased  estimators.  An 
unbiased  estimator  whose  covariance  matrix  is  close  to  the  lower  bound 
will  then  have  an  accuracy  proximate  to  that  of  a  "minimum  variance 
unbiased  estimate."  It  is  shown  in  Appendix  A  that  the  covariance  ma¬ 
trix  of  generalized  least  squares  estimates  is  approximately  equal  to 
the  Cramgr-Rao  lower  bound. 

As  applied  to  the  estimation  problem  of  this  section,  the  Cramgr- 
Rao  Inequality  takes  the  following  form.  Let  the  n  x  1  measurement  vec¬ 
tor  m  have  a  multivariate  normal  distribution  with  mean  vector  p  and  a 
positive  definite  covariance  matrix  Z.  Let  F(p,  S)  ■  0  be  an  r  x  1 
vector  of  constraints  relating  the  n  x  1  measurement  parameter  vector 
p  to  the  p  x  1  estimation  parameter  vector  3.  Assume  that  F(p,  0)  has 
continuous  partial  derivatives.  Suppose  p  £  r  s  n,  and  define 


where  F  and  F0  are  evaluated  at  the  true  values  of  p  and  3. 
u  p 


Let  E"  be  the  covariance  matrix  of  an  arbitrary  unbiased  estimate 

P 

of  6.  Then,  assuming  that  the  regularity  conditions  given  in  Appendix 

B  hold  and  that  the  rank  of  F  is  r  and  the  rank  of  F0  is  p,  we  have 

U  ts 

from  Appendix  B, 

ES  *  Ee  !  <1s(VF*)'1V'1 

Equation  (6)  is  the  Cramdr-Rao  inequality*  for  the  covariance  matrix 

of  unbiased  estimates**  of  £,  and  E.  is  the  Cramdr-Rao  lower  bound. 

P 

The  inequality  E£  >  means  that  the  matrix  E*  -  Eg  is  non¬ 
negative  definite.  This  implies ,  for  example,  that  the  variance  of 
any  unbiased  estimate  of  an  individual  component  parameter  of  B  is 
larger  than  the  corresponding  diagonal  element  of  Eg.  Also,  if  the 
component  estimates  are  approximately  normally  distributed,  the  small¬ 
est  volume  with  a  specified  probability  content  will  be  larger  for  a  co¬ 
variance  matrix  from  Eg  than  from  Eg. 

Equation  (6)  can  be  written  in  alternative  computational  forms  in 
special  cases.  From  Eq.  (6) . 

£8  -  vvv'V'1 

-  <(AFB)*(AFuZ(AFy)*)"1AFB)"1  (7) 

where  A  is  any  nonsingular  r  x  r  matrix.  Thus,  for  example,  in  the 
"just-determined"  case  where  r  ■  p  we  may  set  A  -  Fg^  so  that 

E8  ’  fB1fuEP!<,81>*  <8) 

*Appendix  B  contains  a  proof  of  a  more  general  version  of  the 
Cramdr-Rao  equality  in  which  unbiased  estimates  of  both  p  an  d  B  are 
considered. 

**It  is  shown  in  Appendix  A  that  the  covariance  matrix  of  gener¬ 
alized  least  squares  satlmates  is  approximately  the  Cramdr-Rao  lower 
bound. 
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Equation  (8)  ia  equivalent  to  the  classic  "propagation  of  error"  vari¬ 
ance  relationship  given  in  many  texts  on  applied  analysis.  If  r  -  n 
and  Fy  is  nonsingular,  we  may  set  A  •  in  Eq.  (7)  to  obtain 

£s  -  <»> 

Equations  (7),  (8),  and  (9)  may  be  evaluated  directly  in  terms  of 
F  and  F„  or  in  the  following  manner.  Taking  the  total  differential 

M  P 

of  F(p,  8)  -  0,  we  obtain  the  set  of  simultaneous  equations 

Fydp  +  Fgdfl  -  0  (10) 

Any  permissible  linear  operations  on  Eq.  (10)  to  simplify  it  are  equiva¬ 
lent  to  multiplying  by  a  nonsingular  matrix  A  so  that 

AF^dy  +  AFgd$  -  0  (11) 

Comparing  the  coefficients  in  Eq.  (11)  to  the  expressions  in  Eq.  (7), 
we  see  that  if  Eq.  (11)  is  written  in  an  equivalent  form 

Bdu  +  CdB  -  0 


then 


Eg  -  (C*(BEB*)_1C)_l  (12) 

For  example,  if  Fg  is  nonsingular,  and  Eq.  (10)  is  solved  for  dB  so 
that 


dB  “  Bdy 


then 


BEB 


(13) 


-10- 


Suppoae 
ramete  re  if  necaaaary)  as 


that  F(p,  B)  and  Z  can  be  wtitten  (by  reordaring  the  pa- 


F(p,  B) 


o.  e2>' 

P2(0.  w2.  o.  b2)  J 


(14) 


r  - 


ri  0 


o  z. 


(15) 


(Many  of  the  problems  in  Section  IV  satisfy  Eqs.  (14)  and  (15)  with  &1 
the  emitter  coordinates  and  02  the  aircraft  coordinates.)  Then  Eq. 

(6)  reduces  to 


.“1 


'“i 

0 


F*  A~^F 

F1Bx  1  1BX 

a  -l 

F.a  A/F 


1* 

F  I,F?  0 

iv.fl  lvx 

-1 

'“i  fie2 

0  F2v2S2F2y2 

0  r2.2 

— 

- 

?*  A-  ^7 


-1 


*  .-1-. 

F„„  A/F, 


*  . -1_ 
+  F,_  A/F, 


1B2A1  F1Bx  F1B2A1 'u2  T  r2B2  2  '2B2j 


(16) 


where  A,  -  F,  If*  and  A,  -  Fa  E^.  Inverting  Eq.  (16)  by  the 
for  partitioned  matrices <*>  yields  (if  F^  is  nonsingular) 


rule 


rB2  "  (f2e2A2lr252)  1 
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so  that 


Vl'Il./ul  +  Pl61FlB2£82F162F«i 


(17) 


<r262  <F2h2I2F2ii2>"1f2S2>"1 


(18) 


Equations  (17)  and  (18)  ara  thus  equivalent  to  the  Intuitively 
reasonable  result  that  we  may  estimate  ^  from  the  measurements  on  y2 
alone  and  then  estimate  6^  from  the  measurements  on  and  the  estimate 
of  e2<  However,  this  result  is  not  true  in  general.  For  example,  If 
^  it  overdetermined  (in  which  case  is  not  a  square  matrix) ,  then 

the  measurements  on  y^  can  contribute  to  improved  estimates  of  B2  *nd 
a  smaller  covariance  matrix  than  that  of  Eq.  (18). 

When  n  2  r  >  p  (the  so-called  overdetermined  case) ,  it  is  some¬ 
times  useful  to  introduce  r  -  p  extra  parameters  into  8  (by  defining 
Bp+p  •••»  8r  equal  to  a  subset  of  Up  •••,  yft,  for  example)  so  that 

the  new  F.  is  nonsingular.  This  is  particularly  convenient  if  the 
P  -1 

resultant  matrix  for  FQ  F  can  be  determined  easily  so  that  Eq.  (8), 

p  y 

which  Involves  no  further  matrix  inversions,  can  be  used.  The  new 
parameter  vector  is  B  -  (Bp  •••,  Bp,  Bp+p  •••»  &r)  and  the  Cramdr- 
Rao  lower  bound  for  the  covariance  matrix  for  unbiased  estimates  of 

ft 

(Bp  Bp)  is  the  upper  left-hand  corner  of  the  new  Eg, 

ADDITIVITY  OF  INFORMATION  MATRICES  FOR  INDEPENDENT  MEASUREMENTS 
The  matrix 


«  -  Fg(FpEF*)-1Fg  (19) 

in  Eq.  (6)  is  the  "information"  matrix  for  the  measurement  vector  m 
relative  to  the  independent  parameter  vector  8. 

Suppose  that  m  consists  of  two  stochastically  independent  compo¬ 
nent  vectors  m^  and  m2  with  moan  vectors  yp  y^,  cover lanes  matrices 
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I,,  E„ ,  and  independent  constraint  relations  F, (n1 ,  0)  -  0, 
F2('j2»  8)  *  0.  Then,  for  the  combined  measurement  vector  m, 


£ 


0 


<V  F20> 


F 


P 


Substituting  in  Eq.  (19)  and  expanding. 


*  ■  Fu<Vi<,1>'lris  +  !niT2V2hrl/\i 


■»  4.  +4 

*1  z 


(20) 


so  that  for  independent  measurements  and  constraints,  information  ma- 
trices  are  additive.  Furthermore, 


Ep  -  r1  -  (*1  +  ♦2)_1 


(21) 


zi6  "  zi8  (Eie  +  lj:i0 


(22) 


from  Ref.  2,  p.  29,  where  and  Ejg  are  the  matrices  in  the  Cra* ir~ 
Rao  inequality  estimates  of  0  corresponding  to  the  measurement  vectors 
m^  and  m2>  respectively. 

In  the  computer  simulation  described  in  the  Preface,  Eq.  (21)  pro¬ 
vides  an  efficient  means  for  accumulating  location  accuracy  data  on  in¬ 
dividual  emitters  when  a  number  of  Independent  measurements  are  taken 
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during  the  course  of  the  simulation.  Only  the  distinct  elements  of 
the  information  matrix  for  the  emitter  location  accuracy  need  be  saved 
in  the  computer  for  each  emitter.  As  additional  measurements  are  made , 
e.g.,  bearing  measurements  from  a  single  aircraft ,  the  information  ma¬ 
trix  is  updated  by  means  of  Eq.  '(20).  At  the  end  of  the  simulation 
run,  the  final  information  matrix  for  each  emitter  is  inverted  to  ob¬ 
tain  Ea,  and  then  Eq.  (24)  of  the  next  section  Is  used  to  obtain  the 
p 

location  CEP. 

THE  EMITTER  LOCATION  ERROR  CEP 

If  the  measurements  have  a  multivariate  normal  distribution,  with 
mean  equal  to  the  measured  quantity ,  the  generalized  least  squares  es¬ 
timate  of  0  will  be  approximately  normally  distributed  with  mean  equal 
to  the  true  value  of  0  and  covariance  matrix  given  by  Eq.  (6)  under 
rather  general  assumptions  (see  Appendix  A).  We  shall  assume  that  the 
location  coordinate  estimates  are  approximately  normally  distributed. 

One  measure  of  location  accuracy  is  the  size  of  the  minimum  area 
having  a  specified  probability  of  containing  the  estimate.  From  nor¬ 
mal  distribution  theory,  this  is  the  area  in  an  elliptical  region 
around  the  true  location  with  area 


-2wa  a  (1  -  p2)*  In  (1  -  P)  (23) 

o  yo 

2  2  - 
where  c„  and  ov  are  the  variances  of  the  coordinate  estimates  x  ,  y  ; 

xo  ”0  o’  Jo 

p  is  the  correlation  between  the  estimates ;  and  P  is  the  specified 

probability.  This  measure  of  accuracy  has  disadvantages.  If  the  var¬ 
iances  are  quite  dissimilar  in  magnitude  or  if  the  estimates  are  highly 
correlated,  the  ellipse  is  quite  elongated  so  that  while  the  area  may 
be  small,  the  error  in  one  direction  may  be  quite  large. 

As  a  seemingly  fairer  single  measure  of  accuracy,  the  CEP  of  the 
location  estimate  is  comnonly  used.  The  CEP  can  be  found  from  stan¬ 
dard  tables,  as  a  function  of  0,0,  and  p.^  When  p  -  0,  so  that 

*0  '  o 

the  coordinate  estimates  are  uncorrelated,  an  empirical  approximation 
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J«.  +»j  .,iA 


The  arithmetic  average  of  Eqs.  (28)  and  (29), 


CEP  •  O.75(og  + 

•  0.75 (a2  +  o2 

'  xo  v 


approximates  the  true  CEP  to  within  10  percent. 
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III.  LOCATION  BY  DIRECTION  FINDING  (DF) 

Direction  finding  techniques  for  passive  emitter  location  utilize 
directional  bearing  measurements  to  the  emitter  from  two  or  more  loca¬ 
tions  along  the  reconnaissance  aircraft's  flight  path.  The  emitter 
location  estimate  is  then  taken  to  be  some  measure  of  the  center  of 
the  intersection  points  of  the  bearing  lines  (actually  curves)  formed 
by  the  intersection  of  the  bearing  measurement  planes  and  the  earth's 
surface. 

A  drawback  to  location  by  direction  finding  from  a  single  aircraft 
is  that  for  accurate  locations  the  emitter  must  be  turned  on  suffi¬ 
ciently  long  for  two  or  more  bearing  measurements  to  be  made  some  dis¬ 
tance  apart.  One  way  to  overcome  this  drawback  is  to  take  elevation 
angle  measurements  as  well  as  bearing  measurements.  However,  as  will 
be  shown,  location  estimates  from  a  combination  of  a  single  bearing 
measurement  plus  an  elevation  angle  measurement  are  as  accurate  as  two 
well-spaced  bearing  measurements  only  when  the  ground  range  to  the  emit¬ 
ter  is  of  the  same  order  as  the  aircraft  altitude,  and  the  relative  ac¬ 
curacy  becomes  progressively  poorer  as  the  ground  range  increases. 

It  will  be  shown  that  when  an  emitter  remains  on  sufficiently  long 
for  multiple  bearing  measurements  to  be  made,  the  CEP  decreases  as  the 
inverse  square  root  of  the  number  of  Independent  measurements  taken. 
However,  if  the  measurement  errors  are  correlated,  which  is  likely  when 
the  measurements  are  taken  at  short  Intervals,  there  is  a  limit  to  the 
"equivalent"  number  of  Independent  measurements  that  can  be  obtained 
(see  p.  32)  and  hence  a  lower  limit  to  the  CEP  achievable  by  taking 
multiple  bearing  measurements. 

A  BRIEF  DISCUSSION  OF  THE  MEASUREMENT  ERRORS 

A  comprehensive  error  analysis  for  a  proposed  DF  system  would  at¬ 
tempt  to  account  for  all  sources  of  error  which  contribute  to  the  final 
emitter  location  estimate  error.  The  sources  of  error  depend  upon  the 
detailed  design  characteristics  of  the  system.  A  discussion  of  these 
sources  is  beyond  the  scope  of  this  report. 
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The  errors  may  be  broadly  categorized  as  either  aircraft  naviga¬ 
tion  errors  or  bearing  measurement  error*.  The  aimraft  nayigetior. 
errors  result  in  errors  in  the  estimated  aircraft  location  and  the  es¬ 
timated  direction  of  the  base  line  from  which  the  direction  bearings 
to  the  emitter  are  measured.  The  bearing  measurement  errors  are  com¬ 
binations  of  errors  in  the  sensor  which  measures  the  angle  of  arrival 
of  the  electromagnetic  wave  and  distortions  in  the  electromagnetic 
wav*  due  to  interactions  with  the  aircraft  structure ,  wave  propagation 
anomalies,  etc.  The  latter  errors  may  be  functions  of  frequency,  di¬ 
rection  of  arrival,  elevation  angle,  and  aircraft  attitude. 

One  simplistic  categorization  of  the  errors  is  given  below.  Fixed 
errors  are  those  which  may  be  considered  to  be  constant  throughout  a 
flight.  Random  errors  may  be  correlated  from  bearing  measurement  to 
bearing  measurement,  but  have  mean  zero. 

1.  Navigation  Errors 

a.  Fixed  navigation  position  error 

b.  Fixed  navigation  heading  error 

c.  Random  navigation  position  error 

d.  Random  navigation  heading  error 

II.  Bearing  Measurement  Errors 

a.  Fixed  bearing  error 

b.  Random  bearing  error 

If  the  fixed  errors  can  be  assumed  to  be  random  errors  which  are 
fixed  for  a  given  flight,  then  these  errors  may  be  treated  as  additive 
to  the  random  errors  at  any  given  bearing  measurement  and  are  common 
to  all.  The  net  effect  in  this  case  is  that  a  given  error  type,  bear¬ 
ing  measurement  errors,  for  example,  can  be  treated  as  random  corre¬ 
lated  errors. 

Certain  of  the  errors  may  be  combined  for  error  analysis  purposes. 
The  fixed  navigation  position  errors,  la,  will  simply  add  to  the  emit¬ 
ter  location  errors  if  they  are  constant  in  each  coordinate.  The  fixed 
navigation  heading  error,  lb,  and  the  fixed  bearing  error,  Ila,  may  be 
added  together  as  a  single  fixed  bias.  The  random  navigation  heading 
error,  Ic,  and  the  random  bearing  error,  lib,  may  be  added  together  as 
a  single  random  bearing  error. 


THE  FUNCTIONAL  EQUATIONS  RELATING  THE  PARAMETERS 

The  coordinate  system  used  for  DF  systems  is  generally  some  form 
of  earth  surface  coordinates  and  altitude.  Tor  accurate  location  es¬ 
timation  procedures,  correction  for  earth  curvature  muBt  be  made  when 
the  surface  coordinate  aystem  is  approximated  locally  by  a  rectangular 
coordinate  aystem.  However,  aince  the  earth  central  angle  between  air¬ 
craft  and  emitter  la  quite  small,  the  corrections  are  minor  and  may  be 
considered  equivalent  to  removing  biases  in  the  eatlmatea.  Thus,  for 
location  error  analysis  purposes  it  is  adequate  to  assume  a  local 
rectangular  coordinate  system  tangent  to  the  earth's  surface  in  the 
vicinity  of  the  aircraft  and  the  emitter.  Let  x,  y  be  rectangular  co¬ 
ordinates  in  the  tangent  plane  and  z  be  the  altitude  above  that  plane. 
Denote  the  aircraft  location  by  x,  y,  z  and  the  emitter  location  by 
xg,  yQ,  as  indicated  in  the  following  sketch. 


z 
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Let  9  be  the  directional  bearing  to  the  emitter  from  the  aircraft, 
as  projected  onto  the  x-y  plane.  Since  the  aircraft  and  emitter  are 
not.  In  general,  coaltitude,  the  actual  measured  bearing  angle  is  in 
a  plane  determined  by  the  locations  of  the  aircraft  and  the  emitter 
and  the  axis  of  the  DF  antenna  system  which  we  take  to  be  the  bearing 
base  line  from  which  the  directional  bearings  are  measured  (cf.  the 
preceding  sketch).  With  no  loss  in  generality,  assume  (temporarily) 
that  the  bearing  base  line  is  parallel  to  the  x-y  plane. 

Let  $  be  the  actual  measured  bearing  angle.  Then 

cos  *  "  |vj|R|  (31a) 

where  V  is  a  vector  along  the  bearing  base  line,  R  is  the  vector  from 
the  aircraft  to  the  emitter,  and  VaR  is  the  vector  dot  product. 

Let  i,  j,  and  k  be  unit  vectors  in  the  directions  of  the  orthog¬ 
onal  coordinate  axes.  Then 

R  -  (xq  -  x)l  +  (yQ  -  y) J  +  (zQ  -  z)k 

and 


V*  (R  -  <zQ  -  z) k) 
|V||R-  (,0  -  ,)k| 


<31b) 


since  V  is  perpendicular  to  k.  Thus,  from  Eqs.  (31a)  and  (31b), 


cos  $ 
cos  8 


I*  -  <*o  -  *>*1  JJ-  .2 


■  cos  a 


where  p  ■  |r|,  a  -  z  -  z,  and  a  is  the  elevation  angle  from  the  emit¬ 
ter  to  the  aircraft.  Equation  (32)  is  valid  for  anv  bearing  base  line 
in  the  constant  altitude  plane. 
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Now  assume  that  the  bearing  base  line  Is  parallel  to  the  x-axis. 
We  have 


tan  0  “ 


(33) 


2 

cos  a 


cos  $  ■  cos  8  cos  a 

-  a2  (XQ  -  x>2+  fr0  -  ^ 

P2  <x0  -  x)2  +  (yo  -  y)2  +  (zQ 


(34) 

(35) 


Equations  (33),  (34),  and  (35)  are  the  primary  equations  relating  the 
bearing  angle  <j>  and  the  elevation  angle  a  to  the  aircraft  and  emitter 
locations . 

The  contribution  of  altitude  errors  to  the  total  emitter  location 
error  is  generally  quite  small.  This  can  be  indicated  as  follows: 

From  Eqs.  (34)  and  (35), 


cos  6  -  cos  t 


-sin  0d0  -  coa 


B) 

B) 


-1/2 

-3/2 


da 


-  cos  0  a(p2  -  a2)-1  da 
-2 

d9  ••  -cot  6  ar  da 


(36) 


where  r2  ■  (p2  -  a2)  ■  (xQ  -  x)2  +  (y  -  yQ)2.  From  Eq.  (33), 

2  y0  -  y  1 

sec  9d9  -  - -  dx  -  -  — -  ■  dy 

(xQ  -  x)2  0  Xo  X  0 


y  -  y  x  -  x 

“9  '  -*~T-  d\  '  -V 


dy„ 


■  —  (sin  0  dx  -  cos  6  dy  ) 


(37) 


Substituting  Eq.  (37)  in  Eq.  (36)  and  simplifying. 


~  da  ■  sin  6  cto  -  tan  0  sin  9  dx^  (38) 

r  o  o 

The  bearing  angle  6  is  generally  limited  to  an  Included  angle  of  90 
deg  or  less  centered  on  the  perpendicular  to  the  aircraft  flight  path. 
Taking  0  ■  45  deg  as  a  typical  worst-case  condition  and  Betting  dxQ  -  0 
and  dyQ  ■  0  in  turn,  we  obtain,  from  Eq.  (38),  the  first-order  worst- 
case  effect  of  altitude  errors, 

|dy0l  ■  ldx0l  *  ^  If  dal  < 39> 

Thus,  the  first-order  effect  of  altitude  error  is  directly  proportional 
to  a,  the  difference  in  aircraft  and  emitter  altitudes,  inversely  pro¬ 
portional  to  r,  the  ground  range  between  aircraft  and  emitter,  and 
directly  proportional  to  da,  the  altitude  difference  measurement  error. 
Since  altitude  difference  errors  will  generally  be  hundreds  of  feet  or 
less,  the  contribution  of  altitude  errors  to  the  total  error  will  be 
relatively  small  whenever  the  aircraft- to-emit ter  ground  range  is  much 
greater  than  the  altitude  of  the  aircraft. 

Since  the  primary  use  of  aircraft-to-emitter  altitude  difference 
is  in  finding  the  projected  angle  0  in  the  x-y  plane  from  the  measured 
bearing  angle  and  since  altitude  measurement  errors  are  generally 
of  secondary  importance,  the  altitude  coordinates  will  be  suppressed 
in  the  following  development. 

We  shall  next  derive  analytic  results  for  the  emitter  location 
CEP  in  t'te  special  case  where  errors  in  the  bearing  measurements  from 
the  aircraft  to  the  emitter  dominate  the  aircraft  location  errors. 

Then  the  general  procedure  to  be  used  when  both  bearing  measurement 
accuracy  and  aircraft  location  accuracy  contribute  substantially  to 
the  emitter  location  CEP  will  be  illustrated  by  some  sample  calcula¬ 
tions. 
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EMITTER  LOCATION  CEP  FROM  BEARING  MEASUREMENT  ERRORS  ONLY 

The  emitter  location  CEP  will  depend,  in  general,  upon  the  air¬ 
craft  location-error  and  bearing-measurement-error  covariance  matrices 
and  the  relative  locations  of  the  emitter  and  the  aircraft  when  the 
bearing  measurements  are  made.  In  this  subsection  analytic  results 
will  be  obtained  for  the  contribution  of  independent  bearing  measure¬ 
ment  errors  to  the  emitter  location  CEP,  assuming  that  the  aircraft 
locations  are  known  with  negligible  error. 

Let  N  independent  bearing  measurements  8™,  ***,  0™  be  taken  on 
the  true  bearing  8^,  8N  at  locations  (x^,  y^),  (xN,  yN> 

along  the  aircraft  flight  path.  Let  a ^  be  the  common  bearing  measure¬ 
ment  standard  deviation.  Let  the  bearing  measurements  be  taken  rela¬ 
tive  to  a  bearing  line  parallel  to  the  x-axls  so  that  we  have  the 
situation  depicted  below: 


From  Eq.  (33),  the  functional  relations  between  the  true  bear¬ 
ings  and  the  emitter  location  are 

-  tan~:L(y0  -  y ±)/(*0  -  x±)  i  -  1.  •••»  N  (40) 

so  that,  taking  the  total  differential  and  simplifying, 

yn  -  y,  x  -  x. 

dei - -  dxo - ^  dyo  (41) 

ri  ri 

■  ^-(sin 


8idx0  -  cos  9idy0) 


(42) 


-  *****'*'  ’ 
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2  2  2 
where  r1  -  <xQ  -  x^  +  (yQ  -  y±)  . 

Identifying  0  ■  (^ ,  -  -  * ,  6^)  with  p  and  (x  ,  yQ)  with  0  in  the 

development  of  Section  II,  from  Eq.  (41),  is  the  identity  matrix  and 

-<y0  -  y !>/*£  (x0  -  x^/rj 

Fe  "  :  i  (43) 

~(yo  ~  V^N  <xo  “  *N^rN 

so  that,  from  Eq.  (9),  the  Cramir-Rao  lower  bound  for  the  covariance 

matrix  of  any  unbiased  estimate  of  B  ■  (x  ,  y  )  is 

o  o 

SB  ‘ 

t(y0  -  j,)2/'}  -H\  -  y,X*o  - 

m 

-S(y0  -  y*)^  -  x^/rj  £(x0  -  x^/r* 

(44) 

Two  Bearing  Measurements 

Chooso  the  coordinate  system  so  that  the  x-axis  is  parallel  to  the 
line  segment  joining  the  two  aircraft  positions  and  let 

h  -  y0  -  h  -  y0  -  y2 


Then,  Eq.  (44)  can  be  written  as 
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Since 

[*  c 

Lc  b_ 

the  variances  and  covariance  in  are 

p 

o2  -  h2a^>-1(«ln2  ex  C082  ex  +  sin2  62  cos2  62)  (47) 

o 

o2  -  h2a2D_1(sin4  e±  +  sin4  62>  (48) 

ax  -  h2a2D  1(sin2  0^  cos  0^  +  sin^  0„  cos  02)  (49) 

oJo 

where  D,  after  simplifying,  is 

D  -  sin2  0^^  sin2  «2  sin2  (02  -  0^  (50) 

Substituting  Eqs.  (47)  through  (49)  into  Eqs.  (26a)  and  (26b),  and 
simplifying, 

2  2  2  -1  2  2 
a*  «  h  o*D  [sin  0^  +  sin  02 

+  (sin4  01  +  sin4  02  +  2  sin2  0^  sin2  ©2  cos  2(62  -  6^) )£] 

(51) 

a2  -  h2a2D-1[sin2  02  +  sin2  02 

-  (sin4  0^  +  sin4  02  +  2  sin2  sin2  ©2  cos  2(02  -  0^))^] 

(52) 

Figure  1  contains  CEP  isocontours  as  a  function  of  0  -  0^  and 
A0  ■  02  -  ei  as  obtained  using  og  anc  from  Eqs.  (51)  and  (52)  in 
Eqs .  (24a)  and  (24b). 


i-l 


-  (ab  -  c2)-1 


b  -c 
-c  a 


(46) 
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From  Eq.  (30),  the  CEP  may  be  approximated  as 


CEP  -  0.75 


wherf  from  )’qs.  (47)  and  (48), 


(53) 


°x  +  °y  "  h2ag(sin2  +  »in2  e2)/(ain2  0X  sin2  62  sin2(62  -  6^) 


The  minimum  value  of  the 
@2  ■  120  deg,  where 

min  CEP 


approximate  CEP,  Eq. 

-  /2*  hcf0  (o0  in 

•  0.025  ho0  (o0  in 


(53),  is  at  0^  ■  60  deg, 


radians) 

(54) 

degrees) 

(55) 

Multiple  Bearing  Measurements 


When  the  aircraft  is  flying  a  straight  path  and  a  number  of  bear¬ 
ing  measurements  are  taken  at  nearly  equally  spaced  intervals  along  the 
flight  path,  a  simple  approximation  to  the  location  CEP  can  be  derived. 
For  convenience  in  the  derivation,  assume  a  coordinate  system  as  indi¬ 
cated  below: 


where  the  aircraft  is  flying  in  the  direction  indicated,  h  is  the  per¬ 
pendicular  distance  from  the  emitter  to  the  flight  path,  s  is  negative 


r><u« 
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to  the  left,  and  0^,  0^,  •*•,  0^  are  the  true  directions  from  the  air¬ 
craft  flight  path  to  the  emitter  at  locations  s1 ,  s9,  ••«,  sM  along 
tne  flight  path. 

Identifying  s±  with  -  xq  and  h  with  -  yQ  in  Eq.  (44),  and 
defining  L  -  s^  -  s^,  46  *  6^  -  6^,  and  {1^}  as  the  elements  in  the 
"Information"  matrix  on  the  right-hand  side  of  Eq.  (44),  we  have 

ae  rn  "  E  h2/(h2  +  “2^2 

-  I  /“N  h2/(h2  +  s2)2  ds 

81 

“  ^8h/<h2  +  ®2)  +  tan_1 

“l 

0 

“  2Lh  «in  20  +  6  -  f  ]  N 

®1 

“  2^-  (if  Sin  29jj  -  £  sin  20^^  +  A6) 

“  2Lh  <8in  A0  008  (61  +  V  +  A9)  (56) 

Similarly, 


°02  *22  "  Z  8i2/<hZ  +  8i>2 


“  2lh  ^  “sin  A6  C0B  ^61  +  9n^  + 


(57) 
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ffg  I12  «  “S  h8l/(h2  +  82)2 

=  -  £  /  N  hs/(h2  +  s2)2  ds 

a 

a  JL  Bin  A0  sin  +  6^)  (58) 

From  Eqs.  (56)  through  (58)  and  Eq.  (46)  It  follows  that 


A0  -  sin  46  cos  (6^  +  9^)  2 

N  (40) 2  -  sin2  46  ° 


(59) 


46  +  sin  40  cos  (0^  +  0^)  2 

(40)  -  sin  40 


2Lh  sin  40  sin  (02  +  ©N)  2 

N  (46) 2  -  sin2  49  °9 


(60) 

(61) 


The  rotation  to  a  coordinate  system  in  which  the  estimates  of  the 
emitter  coordinates  are  Independent  is,  from  Eq.  (25),  through  the 
angle 


a 


i  tan-1 


-  4  tan”1 


2  sin  46  sin  (0^  +  8^) 
2  sin  46  cos  (6^  +  6^) 


-  *  +  v 

which  shows  that  the  direction  of  largest  'variation  in  the  joint  distri¬ 
bution  of  estimated  emitter  coordinates  is  in  the  direction  (0,  +  3„)/2 

1  N 

relative  to  the  flight  path. 
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Substitutlng  Eqs.  (59)  through  (61)  into  Eq.  (27)  yields 

2Lh  2  48  +  [4  sin2  A8  cos2  (81  +  8N)  +  4  sin2  48  sin2(81  +  8^) 2 
1  2(46  +  sin  46) (48  -  sin  48)  °8 


.  2UI _ 1 _  2 

H  46-  sin  48  8 


Similarly,  from  Eq.  (26), 


2  2Lh  1  2 

S  N  48  +  sin  48  8 


(62) 


(63) 


Since  L  «  s^  -  s^  «  h  cot  eN  “  h  cot  6^,  we  have 


h  /2  (cot  eN  -  cot 
°Ii  (  68  —  sin  66  J  ® 8 

h  (2(cot  eN  - cot  'A* 

S  ^  l  48  +  sin  48  J  8 

Ik  a  / 48  +  sin  Ae\^ 
og  l 68  -  sin  48 J 


(64) 

(65) 

(66) 


Figure  2  shows  how  the  ratio  o^/Og  in  Eq.  (66)  varies  as  a  func¬ 
tion  of  the  emitter  Included  angle  48  ■  0N  -  8^. 

Figures  3a  and  3b  contain  values  of  1c  =  k(6^,  48)  obtained  using 
Eqs.  (64)  and  (65)  in  the  CEP  approximation,  Eqs.  (24a)  and  (24b)  ,  so 
that 

ho0 

CEP  -  k  — - 

/n 


where  og  is  measured  in  degrees.  Figure  3  is  accurate  enough  for  most 
purposes  when  N  is  greater  than  about  5  and  the  bearing  measurements 
are  more  or  less  uniformly  distributed  along  the  flight  path. 


-  9 ,  (deg) 


-3 


~ 0. (deg) 
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Fig.3b— CEP  isocontours  for  multiple  bearing  measurements 
(small  included  angles) 
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MINIMUM  CEP  FOR  BEARING  MEASUREMENTS  TAKEN  AT  A  CONSTANT  RATE 


When  bearing  measurements  are  taken  at  equally  spaced  points  along 
the  flight  path,  the  CEP  decreases  as  each  new  measurement  is  taken. 
However,  the  minimum  value  of  the  CEP  is  not  zero.  In  practical  terms, 
the  lower  limit  is  reached  when  the  emitter  passes  beyond  the  line  of 
eight  from  the  airborne  platform.  A  general  lower  bound  can  be  ob¬ 
tained  ao  follows.  Since  the  CEP  decreases  aa  A6  increases,  and  ap¬ 
proaches  oL  aa  A6  approachsa  its  maximum  of  ir  radians,  we  may  simply 
evaluate  Eq.  (24a)  at  A8  -  it.  We  thus  obtain 


CEP  %  0.59  /2Lh/N  2ir^ 


*  0.016  vyir  o. 


(c0  in  radians) 


(00  in  degrees) 


where  N '  Is  the  bearing  measurement  density  in  measurements  per  unit 
distance  (in  the  same  units  as  h). 

CORRELATED  BEARING  MEASUREMENTS 

When  bearing  measurements  are  taken  on  an  emitter  with  only  small 
changes  in  the  bearing  angle  between  consecutive  measurements,  the  bear¬ 
ing  measurement  errors  may  be  highly  correlated.  The  simplest  type  of 
correlation  structure  to  treat  la  that  arising  when  the  bearing  measure¬ 
ment  errors  form  a  stationary  Markoff  process.  The  Increase  in  the 
location  estimate  CEP  due  to  correlation  in  the  bearing  measurement 
errors  will  be  illustrated  for  this  case. 

Let  the  correlation  between  successive  measurements  be  equal  to 
p.  From  the  aaeumption  that  tha  bearing  maaeurement  errors  font  a 
stationary  Markoff  process, 


l  ■ 


1  P  ...  P1 
P  1 


e  e 

'  P 


P  . . .  pi 


and 

0 


c;2d-p2) 


i  -p 


iV 


O 


O  X 


w*  -o 


-p 


(67) 
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Asauming  a  linear  flight  path.  have. 


"2  O' 

-h/r‘  8l/rJ 
Lh^n  -n^nj 


!s  ■  <r6rVl 


fi«ui  Eqs ,  (4d)  and  (44), 

(68) 

(69) 


2  2  2 

where  r^  -  h  +  s^,  and  h  and  are  defined  aa  in  Che  sketch  on  p.  26. 
Substituting  Eqs.  (67)  and  (68)  into  Eq.  (69),  performing  the  indicated 
matrix  multiplications,  and  simplifying,  we  find  that  for  large  N, 


ze  “  ^eV’1  °\  +  p)/(1  “  p)  <70) 


Comparing  Eq.  (70)  with  Eq.  (44),  we  aee  that  Figs.  2  and  3  may  be  used 
with  cr0  replaced  by  a0/(l+p)/(l-p)  when  the  bearing  measurement  errors 
form  a  stationary  Markoff  process  with  correlation  between  successive 
measurements  equal  to  p. 

BEARING  MEASUREMENTS  PLUS  ELEVATION  ANGLE  MEASUREMENTS 

If  an  elevation  angle  measurement  is  made  along  with  each  bearing 
measurement,  an  emitter  location  estimate  can  be  made  from  measurements 
taken  from  a  single  aircraft  position.  It  is  shown  in  the  following 
that  when  the  elevation  angle  measurement  accuracy  is  of  the  same  order 
as  the  bearing  measurement  accuracy,  the  combined  bearing  measurement  1 

plus  elevation  angle  measurement  la  roughly  equivalent  to  two  optimally 
spaced  bearing  measurements  when  the  aircraft-to-emitter  ground  range  is  f 

equal  to  the  aircraft-to-emitter  altitude  difference.  When  the  ground 
range  to  the  emitter  is  much  larger  than  the  aircraft  altitude,  eleva¬ 
tion  angle  measurements  must  be  much  more  accurate  than  bearing  mea¬ 
surements  to  provide  equivalent  emitter  location  accuracy.  Thus,  com¬ 
bined  bearing  and  angle  measurements  are  most  useful  against  emitters 
that  are  close  to  the  aircraft  and  when  measurements  can  be  made  at 
only  a  few  positions  along  the  flight  path  (e.g. ,  measurements  on  emit¬ 
ters  with  very  low  power  or  intermittent  operation). 


X-rss^r  TTK' 


ui'f'Vft'/  '---  ••-'••vrw.v-v.v 


Let  a  be  the  aircraft-to-emitter.  altitude  difference,  r  be  the 
aircraft-to-emitter  ground  range,  and  a  be  the  elevation  angle  from 
the  emitter  to  the  aircraft.  Than 


tan  a  »  * 


,  .  (x  -  x  y  -  y 

aec  ado  -  -  f— -y  dxo  +  --j-  dyQ 


2  2  2 
since  r  ■  (x  -  xq)  +  (y  -  yQ)  .  Thus 


hi  -  x  y  -  y 

da  ■  sin  o  cob  a  - = —  dx  + - y-  dy 

V  t  rz  0 


From  Eq.  (41), 


y  -  y0  *  -  xo 

r2-  +  r2,  dyo 

r  r 


Identifying  (6,  a)*  with  u  end  (xo,  yQ)*  with  0  In  Section  II,  we 
have,  from  Eqs.  (71)  and  (72), 


\\  ■  7 


-<y  -  y«> 


(x  -  xo)  sin  a  cos  a  (y  -  yo)sin  a  cos  a 


so  that  for  independent  angle  measurements,  from  Eqs.  (9)  and  (73), 


a~2  0 

■  <V*  ,  \\ 

0  c~z 
a 


(y  -  y0)2  +  c(x  -  x0)2  (c  -  l)(x  -  XQ)(y  -  y0) 

<c  -  l)(x  -  x  ) (y  -  y  )  (x  -  x  )2  +  c(y  -  y  )2 


-4  -2 
r  °9 


(74) 
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2  2  2  2 
where  c  -  sin  o  cos  e*  aa/o  . 

6  01 

Nov  choo se  the  coordinate  system  so  that  x  «  x  and.  therefore, 

2  2  ° 
r  ■  (y  _  y  )  ,  Then  Eq.  (74)  can  be  written  as 


Thus, 


r-1 


[::] 


o  -  ro. 
x  9 

o 


(Mg) 


-2 


o  -  we/^ 

J  n 


ra  /sin  a  cos  a 
a 


r2  +  a2 
- i - °a 


(76) 


where  a^  and  oy  are  squivalc.nt  to  the  crosarange  error  and  range 
error  standard  deviations,  respectively. 

From  Eq.  (30),  the  approximate  emitter  location  CEP  from  the  bear¬ 
ing  measurement  and  elevation  measurement  combination  is  then 


CEP 


-  0.75  (a2  +  a2  )£ 


x  y 

On  *0 


m  0.75  ra. 


1  ♦£♦*)* ;! 


9J 


* 


(77) 


When  aa  ■  a0  ■  cr,  Eq.  (77)  can  be  approximated  as 


CEP  “  0.75  a  r  »  a 


“  0.75  /5  ra  r  «  a 


“  0.75  ao 


a  »  r 


(78) 

(79) 

(80) 
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EMITTER  LOCATION  CEP  CALCULATIONS  FOR  THE  GENERAL  CASE 

Consider  the  following  sketch  depicting  an  emitter  location  (xo> 
y  )  and  five  locations  (x^  y.J ,  •••,  (x5,  y$)  along  the  aircraft 
flight  path: 


(*1»  yx)  Cx2,  y2)  •••  (x5,  y5) 


Assuming  the  bearing  measurement  base  line  is  along  the  x-axis , 
we  have  for  an  aircraft  location  (x,  y) 

tan  6  -  (yQ  -  y)/(xQ  "  *> 

sec2  0d6  -  (x  -  x)_1(dy  -  dy)  -  (x  -  x)"2(y  -  y)(dx  -  dx) 

O  O  o  o  u 

r2d0  -  (xq  -  x) (dyo  -  dy)  -  (yo  -  y) (dxQ  -  dx) 
r2d0  +  (xq  -  x)dy  -  (yo  -  y)dx  -  (xq  -  x)dyQ  -  (yQ  -  y)dxo 

(81) 

where  r2  -  (xq  -  x)2  +  (yQ  -  y)2.  Taking  xq  -  0  and  setting  h  -  yo  -  y± 
(assuming  a  linear  flight  path) ,  Eq.  (81)  may  be  written  for  each  air¬ 
craft  location  as 

hdxQ  +  x^dyQ  "  ~ri^®i  +  +  xi^i  i  "  1*  **•»  5  (82) 

Suppose  that  the  aircraft  location  is  estimated  by  a  navigation 

system  with  independent  errors  in  each  coordinate  with  mean  zero  and 
2 

common  variance  end  that  the  bearing  measurement  plus  navigation 


-38- 


heading  errors  have  mean  zero  and  common  variance  o^.  Setting 

6  “  <V  yo^  and  p  “  ^ei*  X1»  yr  “•»  e5»  x5,  y5),  the  Cramdr-Rao 

lower  bound  for  rh*  I'rtVflHflfti'o  ma  ♦>  «>-(  v  €  ...V  4  - j  ...  •  -  - 

- - - -  ““vaaoBU  «=a  L^uidtCB  Ui  p  18, 

from  Er  j), 


h  ■ 


where,  from  Eq.  (82),  the  iC^  component  of  dF(p,  8)  is 
2 

-r1dei  +  hdx±  +  xidyi  -  hdxQ  -  Xldyo  (84) 

so  that 


Using  equally  spaced  aircraft  locations  from  0.  -  45  deg  to  9C  -  135 


0  *2  .6  .8  1.0 
o  (n  ml) 


Aircraft  location  at3ndard 
deviation  o  (n  ml) 


deg,  a  computer  program  to  evaluate 
a  general  version  of  Eq.  (83)  pro¬ 
duced  the  illustrative  results  in 
the  sketch  at  the  left  for  the 
emitter  location  error  CEP.  The 
sketch  illustrates  the  relative 
contribution  to  the  emitter  loca¬ 
tion  CEP  from  bearing  measurement 
errors  and  aircraft  location  errors. 
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IV.  LOCATION  BY  TIME  OF  ARRIVAL  (TOA) 

Time-of-arrival  techniques  for  estimating  emitter  location  use 
the  arrival  time  of  an  emitter  pulse  at  two  different  aircraft  as  a 
primary  measurement.  The  difference  in  arrival  time  multiplied  by  the 
pulse  propagation  velocity  measures  the  difference  In  range  between 
the  emitter  and  the  two  aircraft.  Curves  of  constant  range  differ¬ 
ences  are  approximately  hyperbolae  on  the  earth's  surface,  and  the  emit¬ 
ter  location  is  estimated  as  the  intersection  of  two  such  curves.  Con¬ 
sequently,  three  or  more  aircraft  must  be  used  for  an  instantaneous 
location  "fix"  on  an  emitter.  Only  two  aircraft  are  required  if  two 
or  more  range  difference  measurements  can  be  made  on  the  same  emitter 
over  some  elapsed  time  interval. 

Two  hyperbolae  can  result  in  as  many  as  four  points  of  Intersection. 
When  there  are  multiple  intersection  points,  additional  information  is 
required  to  eliminate  the  "ghost"  Intersections  from  the  one  represent¬ 
ing  the  true  emitter  location.  The  additional  Information  may  be  pro¬ 
vided  by  TOA  measurements  from  new  aircraft  locations ,  a  crude  directlon- 
of-arrlval  measurement,  or  knowledge  of  the  generel  location  of  the  emit¬ 
ter. 


THE  TOA  HYPERBOLA 

For  illustrative  purposes,  consider  the  two-dimensional  flat-earth 
approximation  with  the  aircraft  and  emitter  taken  to  be  in  the  same 
plane.  Assume  a  coordinate  system  in  which  the  two  aircraft  are  one 
unit  apart  on  the  x-axis,  with  the  origin  at  the  midpoint.  Then  the 
range  difference  to  a  target  at  (x,  y)  is 


d  -  [(x  +  $)2  +  y2]*  -  [<x  -  $)2  +  y2]* 


(85) 


For  fixed  d,  Eq.  (85)  reduces  to  the  branch  of  the  hyperbola, 


2  2 
2L_  _  _X _ 

d2  1  -  d2 


d  £  1 


(86) 
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which  passes  through  the  point  y  ■  0,  x  ■  d/2  and  whose  values  are 
asymptotically  equal  to  the  lines  through  the  origin  with  slopes 

m  ■  ±  (1/d2  -  1)^ 

The  tangent  to  the  hyperbola  bisects  the  angle  formed  by  the  lines 
from  the  emitter  to  the  two  aircraft  locations.  This  can  be  shown  as 
follows : 


Let  P  and  p  be  two  points  on  a  given  hyperbola.  Then  by  the  law 
of  cosines 

d2  -  r2  +  (d^  +  x)2  -  Zrtdj^  +  x)  cos  a 


so  that 


cos  ct  -  (d2  -  r2  -  (d^  +  x)2)/2r(d^  +  x) 


Thus, 


lim  cos  a  ■  -  11m  ^ 
P  -*■  p  P  p 


Similarly, 


lim  cos  g  “  -  lim  ^ 
P  p  P  P 


-Al¬ 


so  that 


11m  a  -  11m  (3 
P  p  P  ■*  p 

shoving  that  tha  tangent  to  tha  hyperbola  bisects  the  angle  formed  by 
the  lines  from  the  emitter  to  the  two  aircraft. 


ERROR  ANALYSIS  EQUATIONS  FOR  THREE-AIRCRAFT  TOA 

For  error  analysis  purposes  and  distances  up  to  a  few  hundred 
miles ,  It  Is  generally  adequate  to  assume  a  flat-earth  coordinate 
system  consisting  of  a  rectangular  coordinate  system  in  a  plane  tan¬ 
gent  to  the  earth  at  a  point  in  the  vicinity  of  the  emitter,  with  alti¬ 
tude  above  the  earth's  surface  as  the  third  rectangular  coordinate. 

The  errors  in  using  such  a  coordinate  system  may  be  considered  as 
biases  to  be  corrected  by  the  location  estimation  procedure. 

Lst  the  aircraft  coordinates  be  (x^  y^ ,  z^)  i  -  1,  2,  3  and  the 
emitter  coordinates  be  (xq,  yQ,  zq).  Let  r^  be  the  distance  from 
location  1  to  location  J ,  and 


*  r ,  -  r. 
io  jo 


(87) 


be  the  difference  between  the  distances  from  the  1th  and  jth  aircraft 


locations  to  the  emitter.  Since  £ 


ij 


dlk  +  .  k  ^  i»  J  »  any  two  of 


^12'  5 ]. 3 •  or  ^31  ma>’  be  taken  to  be  the  TOA  range  difference  parameters 
(TOA  time  difference  multiplied  by  the  speed  of  propagation,  about  1 
ft/nanosecond. ) 

The  equations  relating  the  emitter  and  aircraft  locations  are 


rij  “  (xi  “  xj^2  +  ^yi  “  yj>2  +  "  *j>2  (88) 

Taking  differentials  of  both  sides  of  Eq.  (88)  gives 

drij  ‘  cijx<dxi  -  dxj>  +  cijy(dyi  "  dyj>  +  cijz(dzi  "  dzj)  (89> 
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where 

cljx  ■  <*1  - 
clJj  *  (yl  ■  yl,/rli 

cy*  ■  (,i  -  V/rU 

are  the  direction  coeinea  of  the  line  joining  location  1  to  location  j. 
From  Eqa.  (87)  and  (89), 


(90) 

(91) 

(92) 


"  <ciox  -  ciox)dxo  +  cloxdxl  -  cioxdxi 
+  <cioy  -  cloy)dyo  +  cloydyl  *  cioydyi 
+  (cioz  “  cloz)d*o  +  clozd*l  "  ciozdzi  i  -  2,3  (93) 


Equations  (89)  and  (93)  provide  the  basic  relationships  for  the 
analysis  of  TOA  emitter  location  accuracy.  Equation  (93)  relates  er¬ 
rors  in  the  range  difference  measurements  to  the  errors  in  the  loca- 

A 

tion  coordinate  vector  8  -  (xo>  yQ,  zo>  x^,  y^,  z3),  whereas 

Eq.  (89)  relates  errors  in  range  measurements  between  aircraft  to  er- 

A 

rors  in  8  .  Similar  equations  hold  if  ground-station-to-aircraft 
range  measurements  are  included. 

In  the  following,  we  shall  first  obtain  analytical  results  for 
the  emitter  location  CEP  for  cases  in  which  certain  of  the  measurement 
error  accuracies  dominate  the  remainder.  Then,  sample  calculations 
will  be  presented  illustrating  the  general  procedure  for  obtaining  the 
emitter  location  CEP  when  both  the  aircraft  location  accuracies  and 
the  aircraft-to-emitter  TOA  measurement  accuracies  are  important. 


-43- 


EMITTER-TO- AIRCRAFT  RANGE  DIFFERENCE  ERRORS  ONLY 

Suppose  Che  three  aircraft  locations  and  the  emitter  altitude  are 
estimated  with  negligible  error. 

Let 


51  “  612  "  rlo  -  r2o 

(94) 

62  "  513  "  rlo  -  r3o 

(95) 

2  2 

be  measured  with  measurement  error  variances  c^,  o^  and  correlation 
coefficient  p. 

Assume  also  that  the  aircraft  and  the  emitter  are  in  the  same 
geometric  plane  (or  that  the  altitude  differences  are  inconsequential). 
Then  Eq.  (93)  can  be  written  as 

d61  ■  (cos  y2  “  cos  +  (sin  y2  "  sin  Yl^dyo  (96> 

d62  -  (cos  y3  ~  cos  Y]^)^  +  (sin  Y3  "  sin  Y^)^yo  (97) 

where  is  the  angle  between  the  line  joining  the  lt*'  aircraft  and 
the  emitter  and  the  x-axis. 

Now  choose  the  x-axis  to  bisect  the  angles  between  the  lines 
joining  two  of  the  aircraft  and  the  emitter  as  indicated  in  the  sketch 
below: 


Then,  from  Eqs.  (96)  and  (97), 


do^  -  [cos  (a  +  0/2)  -  coa  0/2]dxQ  +  {sin  (a  +  0/2)  -  sin  0/2]dyQ 
-  -2  sin  — sin  -y  dxQ  +  2  sin  cos  °  -y  dyQ  (98a) 

d62  -  {cos  (  -  0/2)  -  cos  6/2]dxo  +  [sin  (  -  0/2)  -  sin  (0/2)]dyo 
*  -2  sin  |-  dyQ  (98b) 

Solving  Eqs.  (98a)  and  (98b)  for  dxQ  and  dyQ, 

dxQ  »  -  esc  n  ■  ^csc  Tf  d51  +  esc  cos  -~£~  dd^ 

dyQ  »  -  esc  ddj  (99b) 

From  Eqs.  (99a)  and  (99b),  after  some  reduction, 


+  cr2  -  $  esc2  -2-—^-  (oj  esc2  |  +  c2  esc2  | 

+  2pc102  CSC  esc  cos  ^ 
*  d’"1^o2  sin2  -|  +  a2  sin2  j  +  2p0jO2  Bin  f  Bil1  2^  C0B  "Y 


where 

D  -  4  sin2  |  sin2  j  sin2 
*  J[sin  a  +  sin  0  -  sin(a  +  0)]2 


(101) 


Substituting  half-angle  formulas  in  Eq.  (100)  and  simplify lug  yields 


2  2 

,  ,  (l-cosoja,  +  (l-cosB)o,  +  pa.c, (cosa  +  cosB  -  1  -  cos(ot4-S)) 

<T  +  m  2 - ± - - - — - - - 


x  y 

o  7o 


(sina  +  sinB  -  ain(o+B)) 


From  Eq.  (96) , 


61  "  rlo  "  r2o  "  tloVlo  "  fc2ov2o 


52  *  rlo  “  r2o  “  CloVlo  ■  t2ov2o 


where  tio  ia  the  propagation  time  and  the  average  propagation  ve¬ 
locity  (about  1  ft/nanoBecond)  from  the  emitter  to  the  1th  aircraft. 

Assuming  that  the  variances  and  covariances  of  the  errors  in  the  mea- 

2 

surements  of  t^v.  are  equal  with  common  variance  and  common  cor¬ 
relation  pd  (which  ia  reasonable  if  either  the  velocity  estimate  errors 
are  negligible  compared  to  the  time-o£-arrival  errore  or  if  the  emitter 
distances  from  each  aircraft  are  comparable) ,  the  covariance  matrix  of 
the  measurement  errors  of  6^  and  6^ 


.pclff2 


2  1 


1  2 


where  c2QA  -  (1  -  pd)cd  is  the  "affective"  range  difference  error  vari¬ 
ance. 

Substituting  Eq.  (103)  in  Eq.  (102)  gives 


_2  .  .2  3  -  cos  a  -  cos  8  -  cos  (a  +  8)  „_2 

o  +  o  - - e 1 2c_ 

o  yo  (sin  a  +  sin  8  +  sin  (a  +  $))* 

3  -  cos  -  cos  n2  “  cos  n3  2 

"  “  ~~~  ~~~  7F~  2°toa 

(sin  +  sin  n2  +  sin  hj) 


(104) 
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where  i  0,  +  n2  +  ^3  “  360  deg,  and  n^»  *  ^3  are  tl'e  angles  be¬ 

tween  the  lines  from  the  emitter  to  the  aircraft,  measured  in  a  clock¬ 
wise  direction. 

Equation  (104)  is  a  minimum  when  n-j_  ■  n2  ■  Dj  *  120  deg,  which  im¬ 
plies  that  the  emitter  location  errors  tend  to  be  smallest  near  the 
center  of  the  aircraft  triangle. 

It  may  be  noted  that  the  emitter  location  CEP  is  proportional  to 
aT0A  an<*  ^dependent  °-  scale,  !•«•»  Eq.  (104)  is  a  function  only  of 
the  angles  between  the  aircraft  and  the  emitter. 

Figure  4  contains  CEP  isocontours  as  a  function  of  a  and  B  (as  S)  , 
the  two  smallest  angles  from  n2>  end  q^.  This  figure  is  thus  valid 
for  all  aircraft  triangles.  Figures  S,  6,  and  7  contain  CEP  isocon¬ 
tours  for  three  isosceles  aircraft  triangles,  as  obtained  from  a  more 
accurate  CEP  calculation  than  the  approximation  of  Section  II.  The 
figures  are  not  valid  for  emitter  locations  naar  the  extensions  of 
the  aircraft  triangle  base  legs,  where  certain  assumptions  made  in 
the  analysis  are  not  met.  In  fact,  Eq.  (104)  implies  an  infinite 
CEP  along  the  base  leg  extensions— which  is  not  true  for  the  location 
estimate  that  is  the  intersection  of  the  range  difference  hyperbolae. 

From  the  figures,  at  a  given  distance  from  the  centroid  of  the  iso¬ 
sceles  triangles,  the  CEP  is  a  minimum  on  the  perpendicular  bisector  of 
the  unit  base  leg.  Along  the  perpendicular  bisector,  ■  n2  ■  n  and 
■  2tr  -  n,  so  that  Eq.  (104)  becomes 


2.2  3-2  cos  n  -  cos  2n  „^2 

°x  +  ay  - - - TX  TOA 

0  Jo  (2  sin  n  -  sin  2n) 


<2nZ  +  <2n)2)/  j^“  - 


aT0A  n  <  <  1 


“  °toa 


*  96  r4»to* 


R  >  >  1 


3 (deg) 
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where  R  Is  the  distance  from  the  emitter  to  the  unit  base  leg.  Thus, 
CEP  *  0.75  i^96  R2o. 


TOA 


«  7  R2o 


TOA 


R  >  1 


2  2 

A  similar  reduction  for  ox  and  oy  for  locations  along  the  per¬ 
pendicular  bisector  yields  0  ° 


«*  *  4  A  8  °10A  R  ~  1 

o 


°y  “  71  80  TOA 

J  A 


R  »  1 


and  the  correlation  term,  by  symmetry,  is  aero.  Considering  the  y-axis 
to  be  "downrange"  and  the  x-axis  to  be  "crossrange," 


°DR/aCR  "  4  ^  R 


■  7R 


R  »  1 


so  that  the  downrange  errors  are  much  larger  than  the  crossrange.  This 
is  also  true  for  locations  off  the  perpendicular  bisector. 

AIRCRAFT  LOCATION  ERRORS  ONLY— I 

Suppose  that  the  aircraft  locations  are  estimated  from  Independent 

navigation  systems  and  that  the  TOA  errors  are  negligible  relative  to 

the  aircraft  location  errors.  Suppose  that  the  aircraft  and  emitter 

altitude  estimate  errors  are  also  negligible.  Assume  that  the  aircraft 

and  emitter  are  In  the  same  plane  (or  that  the  altitude  differences  are 

inconsequential).  Let  the  aircraft  locations  be  (x^,  y^)  1  ■  1,  2,  3 

and  the  emitter  location  be  (xq,  yQ).  Suppose  also  that  the  aircraft 

navigation  system  position  errors  are  independent  in  each  coordinate 

2 

with  common  variance  o  . 
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Slnce  and  15 ^3  are  assumed  known,  from  Eq.  (93)  with  d6j^  “ 
dfi .  m  n  u«  cbtnin 

u 

(cos  y2  •  c°8  Yj_)d*e  +  <®ln  Y2  “  *ln  Yj)dy0 

»  cos  Y2dx2  "  601  Yldxl  +  *in  Y2^y2  "  *ln  Yld5?l 
(cos  Y3  “  cos  Y1)dxo  +  (sin  y2  “  •*»  Y1)dy{) 

■  cos  YjdXj  -  cos  Yjdx^  +  sin  Y^dy^  ~  ®in  Yjdy^  (105b) 

*  / 

Taking  the  measurement  parameter  vector  as  u  *  (x^,  y^,  X2,  y2»  x-j» 
y^),  from  Eq.  (105), 
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independent  ot  the  emit  Lei:- tc  aircraft  range  difference  errors,  from 
Eq.  (104), 


3  -  COS  n,  -  COB  Ho  “  COS  Ho  o  9 

- — - 1 - r  2<ct  +  °toA) 

(sin  +  sin  hj  +  8^n  hg) 


(108) 


for  the  combined  effects  of  the  measurement  errors.  Thus,  Figs.  4 

. 

through  7  may  also  be  used  for  the  combined  errors  by  replacing  oTn, 

2  9  A  1UA 

by  (a  +  <jtoa)^' 

AIRCRAFT  LOCATION  ERRORS  ONLY— II 

As  another  example  in  which  analytic  results  can  be  obtained,  con¬ 
sider  a  TOA  system  in  which  the  aircraft  locations  are  estimated  from 
DME  measurements  between  aircraft  and  aircraft  position  measurements 
from  independent  navigation  systems. 

Suppose  that  the  TOA  measurement  errors,  the  aircraft  and  emitter 
altitude  measurements  errors,  and  the  DME  errors  are  negligible  rela¬ 
tive  to  the  aircraft  position  errors.  Assume  also  that  the  aircraft 
and  emitter  are  in  approximately  the  same  plane  and  denote  the  air¬ 
craft  positions  by  (x^,  i  ■  1,  2,  3  and  the  emitter  position  by 

(x  ,  y  ).  Suppose  also  that  the  aircraft  position  errors  are  indepen- 
0  O  | 

dent  in  each  coordinate  and  from  aircraft  to  aircraft,  and  let  the 

2 

errors  have  common  variance  a  . 

For  convenience,  let  the  coordinate  system  have  its  origin  at  the 
centroid  of  the  aircraft  triangle  so  that  x^  +  Xj  +  x^  *  0  and  y^  +  + 

y3  -  0.  Let  be  the  angle  between  the  x-axis  and  the  line  joining 

(x^,  y^)  and  (xq,  yQ)  and  $  the  angle  between  the  x-axis  and  the  line 
from  the  origin  to  (xQ,  yQ).  Then 


xi  -  xq  -  ri(j  cos  +  A^)  i  -  1,  2,  3  (109a) 


yi  - 


r^o  sin  ($  +  A$j,)  i  ■  1,  2,  3 


(109b) 
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where  4$^  *  $  -  ^  Is  known. 

From  Eq.  (109), 

dxA  »  dxQ  -  rlo  sin  ($  -  Ad> d$ 

-  dxo  -  <yi  -  yo)  d*  i  -  1.  2,  3  (110a) 

dyi  “  dyo  +  rio  COfl  ^  '  A*i>  d* 

-  dyo  +  (x±  -  xo)  d*  i  -  1,  2,  3  (110b) 

& 

Setting  the  measurement  parameter  vector  p  ■  (x^»  Xji  x^,  y^,  y^, 
y3)  and  the  estimation  parameter  vector  0  «■  (xQ,  y  f  $),  from  Eqa.  (110) 
and  (12) , 

Eg  -  (C^C)"1  a2  (111) 

where 

0  y0-?1 

0  V?2 

0  Vy3 
1  “l"*. 

1  Vxo 

1  x3-xo 

Thus, 
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where 


a  - 


2  (y0  -  m  3yc 


i-i 


b  “  j.  <xl  -  xo>  “  ^0 

i-1 


c  -  j  [<xi  "  xo)  +  <yi  -  yQ)  ] 

i-1 


i-1 


<xi  +  #  +  +  yo> 


-  3(r2  +  R2) 

2 

r  1b  the  average  squared  range  between  the  aircraft  locations  and  the 
2 

centroid,  and  R  is  the  squared  range  from  the  emitter  to  the  centroid. 

Now  choose  the  x-axls  to  pass  through  the  emitter  location  so  that 
yQ  -  0.  The  error  in  the  x-directlon  may  then  be  considered  to  be  the 
range  error  and  the  error  in  the  v  direction  the  crossrange  error. 

Then 


and 


* 


C  C  -  3 


1 

0 

.0 


0 

1 

R 


R2+rl| 


0 


0 


R2+r2  -R 


-R 


(112) 
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! 

j 

') 

t 

i 

! 


J 


e  »•  _ 
"JO 


Substituting  Eq.  (112)  in  Eq.  (Ill), 


a l  -  (1  +  R2/r2)o2/3 


(11 2a) 

(113b) 


so  that 


CEP  «  0.75  (a2  +  o2  )£ 
o  yo 

«  0.43  (2  +  R2/r2)&  a  (114) 

When  the  aircraft  triangle  is  equilateral  with  side  .length  A,  the 

_  2  2 

distance  from  the  aircraft  to  the  centroid  is  A// 3,  so  that  r  *  A  /3 
and  Eq.  (114)  becomes 

CEP  ■  0.43  (2  +  3(R/A)2)^  a 

«  0.75  (R/A)  o  R  »  A 

EMITTER  LOCATION  CEP  CALCULATION  FOR  THE  GENERAL  CASE 

In  the  preceding,  we  have  derived  analytic  results  for  the  emitter 
location  estimate  CEP  for  the  TOA  systems  for  cases  in  which  certain  of 
the  measurement  errors  dominated  the  remainder.  The  following  example 
illustrates  the  general  procedure  for  determining  the  emitter  location 
CEP  which  combines  the  contributions  from  all  the  measurement  errors. 
Suppose  that  the  aircraft  locations  are  determined  by  combining  air¬ 
craft  altimeter  measurements  with  DME  range  measurements  between  air¬ 
craft  and  between  the  aircraft  and  two  ground  stations  (as  in  one  mode 
of  the  Air  Force's  QRC-334  system).  Suppose  also  that  the  emitter  al¬ 
titude  is  estimated  separately  by  an  independent  measurement  with  known 
standard  deviation.  Denote  the  (known)  ground  station  locations  by 
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\ 

(x^,  y^,  z^)  and  (x5,  y5,  Zg).  The  parameter  vector  to  be  estimated 
can  then  be  taken  as 

p  -  (*0*  yQ.  zQ»  xi»  yx»  z2*  x2»  y2>z2*  x3*  y3’  z3^ 
and  the  measurement  parameter  vector  as 


* 

P 


613'  r14*  r  15'  r24*  r25'  r34'  r35*  r12'  r13’  r23*  z0* 
From  Eq.  (6),  the  Cramdr-Rao  lower  bound  is 


Z 


B 


F^Z'^)'^)”1 


(115) 


where,  from  Eqs.  (89)  and  (93),  F  -  X,  and  F.  is  the  matrix  on  p.  58. 

M  P 

Suppose  that  the  standard  deviation  of  the  DME  range  measurements 
Is  10  ft,  that  of  the  TOA  range  differences  is  20  ft  with  a  correlation 
of  1/2  (l.e. ,  assuming  independent  time-of-arrival  measurement  errors 
at  the  three  aircraft),  that  of  the  aircraft  altitudes  is  40  ft,  and 
that  of  the  emitter  altitude,  100  ft.  Assuming  that  all  measurements 
except  TOA  range  differences  are  uncorrelated,  the  non-zero  elements 
of  the  measurement  covariance  matrix  £  are  Z^  -  Z22  ■  400,  Z^2  -  Z21  - 
200,  Z33  -  ...  -  EU|11  -  100,  Z12>12  -  10,000,  Z13>13  -  Zw>14  -  E15>15 
1600. 


The  following  figure  illustrates  the  output  of  a  computer  program 
which  calculates  Z^  from  Eq.  (115)  and  then  the  approximate  CEP  of  the 
emitter  location  estimate  from  Eqs.  (26),  (27),  and  (24): 

Numbers  in  a  given  column  are  approximate  emitter  location  CEPa  for 
emitters  at  the  points  indicated  by  the  blocked  number  in  the  column 
when: 

Ground  stations  GS1  and  GS2  only  are  used  and 


a.  Aircraft-to-aircraft  distance  measurements  are  omitted 

b.  All  measurements  are  Included 

c.  There  are  TOA  measurement  errors  only— known  aircraft 
locations 
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Emitter  CEP  (ft) 


Ground  stations  GS3  and  GS4  only  are  used  and 

c.  There  are  TOA  measurement  errors  only — known  aircraft 
locations 

d.  All  measurements  are  included 

e.  Aircraft-to-aircraft  distance  measurements  are  omitted 

As  indicated  in  the  figure  above,  the  CEP  depends  on  the  relative  lo¬ 
cations  of  the  ground  stations  and  the  aircraft,  so  different  results 
obtain  for  each  configuration  of  ground  stations  and  aircraft  triangle 
considered. 


Preceding  page  blank 
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Ap pend lx  A 

GENERALIZED  LEAST  SQUARES  ESTIMATION 


Aa  described  In  Section  II,  the  general  estimation  problem  In  lo¬ 
cation  finding  Involves  a  set  of  measurements  m  ■  (m, ,  • • • ,  m  )  on  the 

e  ^  ^ 

observed  parameter  vector  y  -  •••,  yR)  and  a  set  of  functional 

relationships  between  y  and  the  (possibly)  unobserved  parameter  vector 

_  * 

3  ■  (Bp  *  Bp)  » 

f(y»  3)  -  (f^y.  6),  fr(y,  3))*  ■  o  (116) 


F(y ,  8)  may  Include  relationships  on  y  or  8  alone,  i.e.,  certain  of 
the  f^(y,  B)  may  be  Independent  of  y  or  of  3.  The  general  estimation 
problem  is  to  estimate  y  and  6  subject  to  the  constraints,  Eq.  (116). 

A 

Generalized  least  squares  estimates  of  y  and  3  are  the  values  y 

A 

and  3  which  minimize  the  quadratic  form 

Q  -  (m  -  y)*Z-1(m  -  y)  (117) 


subject  to  the  constraints  F(y ,  8)  ■  0.  In  general,  when  F(y,  8)  is 

A  * 

nonlinear,  the  determination  of  y  and  6  require  numerical  methods.  One 
widely  used  iterative  procedure  is  an  adaptation  of  Newton's  method, 
sometimes  called  the  method  of  linearization. 

Let  y°,  6°  be  initial  guesses  of  y  and  6.  Approximating  the  iC^ 
component  of  F(u,  8)  by  the  first-order  terms  in  the  Taylor  series 
expansion  around  y°,  8°, 


• s)  ■  £i(,,0>  e0>  +  I  -  "j0) 


3f  j 


Z  IT 


k-1 


k'y  -  y  ,  8  -  8 


0(8j  -  B°)  (118) 


I 


S3 
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Deflning  F  as  the  r  x  n  matrix 


p  l  3p  I  o  -O 


and  Fg  as  the  r  x  p  matrix 


"  Kl  -  e  -  e° 


Eq.  (118)  becomes,  in  matrix  form, 


F(p,  8)  -  F(p°,  8°)  +  Fy(p  -  p°)  +  Fg(B  -  8°) 


New  approximations  to  p  and  8  are  obtained  by  first  imposing  the  con¬ 
straint  equation  F(p,  8)  *  0  on  Eq.  (119)  to  obtain 


0  -  F(p  ,  8°)  +  Fy(p  -  p°)  +  Fe(B  -  8°) 


and  then  minimizing  the  quadratic  form,  Eq.  (117),  subject  to  the  lin¬ 
ear  constraints  of  Eq.  (120).  The  latter  problem  is  a  generalization 
of  the  usual  linear  least  squares  estimation  problem.  A  quite  general 
solution  is  given  in  Ref.  4.  When  F  and  F.  have  full  rank,  the  new 

A  A  M  P 

approximations  to  p  and  8  are 


(121a) 


A  more  standard  notation  for  Fy  and  Fg  would  be  F  (p°,  8°)  and 
Fg(p°,  8°)»  Indicating  the  values  of  the  arguments  at  which  the  func¬ 
tions  are  evaluated.  However,  this  is  cumbersome.  The  values  of  p 
and  8  at  which  F  and  Fg  are  evaluated  should  be  clear  from  the  con¬ 
text.  M 
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'  '*r^r*#f&»*i*a<m**~***> 


61  -  -(F*(F  EF*)-1F  )”XF*(F  EF*)_1(F  m  -  v  )  (121b) 

p  y  p  p  B  u  u  v  o 

where  vq  -  f^u°  +  Fg00  -  F(p°,  3°). 

Approximate  covariance  matrices  and  cross-covariance  matrices  of 
p^  and  are 

V  -  E  -  EP?<VpJrl11  " 

(122a) 

V  ■  <rs<Vf*)’1V"1  <222b> 

v»>  ■  <122o) 

If  the  Initial  guesses  p°  and  0°  are  close  to  p  and  0,  minimizing 
the  quadratic  form  of  Eq.  (117)  subject  to  the  linear  constraints  of 
Eq.  (120)  will  be  approximately  equivalent  to  minimizing  the  quadratic 
form  subject  to  F(p,  0)  *  0.  The  result  is  a  new  approximation  p1, 

A  A 

0  to  p,  0.  The  entire  procedure  may  be  iterated  until  convergence 
of  the  sequence  (p*,  0*)  is  indicated.  Conditions  under  which  the 
sequence  converges  to  p,  0  are  contained  in  Ref.  5,  for  example. 

Comparing  Eq.  (122)  with  Eq.  (125)  of  Appendix  B,  we  see  that  the 
only  difference  is  the  value  of  (p,  0)  at  which  and  Fg  are  evaluated. 
When  the  covariance  matrices  of  p  and  0  are  small  so  that  p  and  0  are 
close  to  p  and  0  with  high  probability,  the  approximate  covariance  and 
cross-covariance  matrices  of  Eq.  (122)  will  be  close  to  those  of  the 
Cramdr-Rao  lower  bound  for  unbiased  estimates  of  p  and  0. 
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Appendix  B 

THE  CRAMER- RAO  INEQUALITY 


Let  p(X,  6)  be  the  probebllity  density  function  of  e  vector  rsndou 
variable  X  with  «n  unknown  parameter  vector  6  -  (6^,  . ...  0^)  .  Define 


,  .  _  Ln_E. 

?ij  301  36^ 


♦u  ■  E(pl]> 


The  matrix 

*  -  (♦ij)  U  i  -  1.  ....  k 
is  called  the  information  matrix  for  p(X,  0). 

4* 

Theorem  1.  the  Cran6r-Rao  inequality  for  unbiased  estimators. 
Let  SjCX),  e2(X),  ....  at(X)  be  t  etatietioe  aw oh  that 


X[«£<X)]  ■  •••»  ®jj) 

i»  J  “  1»  •••»  b 

EI(#1  “  gl5(*j  “  “  vij 

Suppose  Sgj/SOj  exists  and  define  the  matrices 


v  *  CVtj)  l,j  -  1,  ....  t 

&  ■  Og^/SOj)  1  ■  1>  ....  t  j  ■  1,  ....  k 

Reference  2,  p.  265. 
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If  differentiation  may  be  passed  under  the  integration  sign  so 

that 


8l(X)p(X,  0)dX 


f  ’i(lt>  a  « 


then 


V  -  A*" A*  i  0  (123) 

i.e.t  the  matrix  V  -  A*  A*  is  nonnegative  definite 

In  Section  II,  the  logarithm  of  the  multlverlate  normal  density 
function  for  the  measurement  vector  m  is  (aside  from  a  term  independent 
of  p)  proportional  to 


“  $  (m  -  w)V1(m  -  u) 


(124) 


where  the  n  x  1  measurement  parameter  vector  y  Is  related  to  the  pxl 
independent  parameter  vector  0  by  the  r  x  1  vector  of  constraints 

F(V,  0)  -  0 


Let  w°,  8°  satisfy  F(p°,  8°)  -  0  and  define 


F 

V  3p 


V°»  0  ■  0° 


(r  x  n) 


3F 

30  U  ■  V°,  0  ■  0° 


(r  x  p) 


*The  matrix  •”  is  any  generalised  Inverse  of  4, 
satisfying  «  ».  If  *  is  positive  definite , 


1  #  e  •  t 
,  *-l. 


any  matrix 


-66- 


Assuming  that  F  has  continuous  partial  derivatives  at  (m  ,  &  ) » 

F  has  full  rank  r,  F0  has  full  rank  p,  n  as  r  i  p,  and  that  the  rsn..|r-- 

U  P 

ir.cr.ts  above  cor  the  Cramdr-Rao  inequality  hold,  we  obtain  the  following 
result. 

Theorem  2.  With  the  above  assumptions,  the  dispersion  matrix  E*  g 
of  an  unbiased  estimate  y.,  $  of  v° t  e°  satisfies  the  matrix  inequality 
(A  >  B  means  that  A  -  B  is  positive  definite) , 


r  a  a 

u.fl 


y  ye 

\  EuB 

m 

i 

fa  *  Pa 

r6y  6 

Eeu  £e 

—  m 

y,e 


(125a) 


where 


\  - 1  -  U25b) 


h  * 


£„a 


Eeu  “  £ws 


(125c) 

(125d) 

(125e) 


Proof.  First,  suppose  that  £  is  the  identity.  The  total  differ¬ 
ential  of  F(y,  6)  “  0  at  u°,  g°  is 


F  dp  +  F.dg  ■  0 
y  p 


(126) 


Let  the  columns  of  the  n  x  n  -  r  matrix  U  and  the  columns  of  the 
n  x  r  matrix  W  be  orthonormal,  and  let  the  columns  of  U  span  the  null 


The  proof  given  here  uses  properties  of  the  pseudo-inverse  of  a 
matrix.  The  pseudo-inverse  of  a  matrix  A  is  the  unique  matrix  A*  sat¬ 
isfying  (a)  AA+A  -  A,  (b)  A+AA+  -  A+,  (c)  (AA+)*  -  AA+,  (d)  (A+A)*  - 
A+A.  The  properties  of  the  pseudo-inverse  used  here  are  contained  in 
Ref.  A,  for  example. 


-67- 


space  of  Fy  and  the  columns  of  W  span  the  space  orthogonal  to  the  null 

«.  £  v 

“r““w  Ay‘ 

Define  6  and  y  by  the  Invertible  transformation 

U  -  W£  +  Uy  +  p°  (127) 


and  H(£,6)  by 


H(6,  8)  -  F(W6  +  Uy  +  y°,  6)  (128) 

A  it  it 

where  0  ■  (8  ,  -y  )  .  Then  H(6,  0)  has  continuous  partial  derivatives , 

A  I  A 

H(0,  eQ)  -  0,  where  8q  »  ($°  ,0)  »  and  ||  L  Q  _  0  -  FyW  is  non- 

x  O 

singular.  By  the  implicit  function  theorem,  H(6,  9)  «  0  has  a  unique 

solution  6  -  h(9),  i.e.,  H(h(0) ,  6)  -  0,  in  a  neighborhood  of  8  .  Thus 

o 

0  is  a  vector  parameter  indexing  the  normal  probability  density  n(m,  p) 
in  that  neighborhood. 

From  Eq.  (127), 


dp  ■  Wd6  +  Udy  (129) 

Substituting  Eq.  (129)  into  Eq.  (126),  we  obtain 

F  (WdS  +  UdY)  +  Fad8  -  0 

M  D 

F  WdS  +  0  +  F0d8  -  0 
so  that,  since  F^W  is  nonsingular, 

d«  -  -  (FuW)_1FBde  (130) 


t 


Reference  6,  p.  138. 
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Substituting  Eq.  (130)  in  Eq.  (129), 


dp  -  -  W(FuW)“1Fgde  +  Udy 

-  -  wwVFgde  +  udY 


or,  In  matrix  fora. 


-  Pj^Fgde  +  tidy 


dp  -  -  (fJfb<  U)d9 


The  Information  matrix  «  la 


-'(siM*) 


where 


±  .  f-L.  _a_Y 

36  l  3B1  .  30  J  * 


We  have 


_  3  In  p  ,3  * 

36  t  30<®  “  W)  (o  - 


P) 


■  *  fe  -  V)*(m  -  p) 

-  (r*r6,  U>*<»  -  M) 


from  Eq.  (132).  Thus, 


(131) 


(132) 


(.133) 


(134) 


-  si  [hH*  ■  ft*  u> 


-  (F+F&.  U>*(F+Fgt  U) 


<<v  *l(<v u) 


’^VX'b  (P>3)*U 

u*f+fb  u*u 


<<VX?s  0 


alnce  U  U  -  I  and 


u*f+  -  oVVy 


Since  the  right-hand  aide  of  Eq.  (135)  la  a  conatant,  It  la  equal 
to  ita  expected  value  and  tharefora  to  *.  It  la  eaay  to  verify  that 


0 


(136a) 


To  apply  the  Cramdr-Rao  Inequality  to  aatlaates  of  u  and  g,  aet 


G(8) 


Then,  In  Theorem  1, 


A  36 

A  m  — 

ae 


Hr  .Hr 


from  Eq.  (132).  From  Eqs.  (136a)  and  (136b), 


* 

-F+F„  -U 
li  6 

1 

O 

ca 

1 

A$  A  ■ 

I  0 

0  1 

L  J 

-u  o 

m  m 

where 


Z 


8 


Performing  the  Indicated  matrix  multiplications  in  Eq.  (137) 


+  * 
A*  A 


<vs(FX>‘  *  w‘  -Krsh 


where,  from  the  definition  of  U,  *U*  ■  I  -  F+F  (see  Ref.  4). 


(136b 


(137) 


ve  have 


(138) 
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where, 

h  ■  <»„*>  V<F„*>V' 

•  <FJv>+*<v)+V+ 

■  XvO'V'1  (143) 

AUuV  -  A(I  -  <PiiA)+(FpA))A* 

-  A(I  -  <F/>*<F/<FyA)VFyA)A* 

-  *d  -  *x<v»>'fy»* 

■  1  -  tP*<Pl,EF‘>"1»w£  (1*4) 

and 

*(FV*)+  . 

"  XVO"1  (1*5) 

Substituting  Eqs.  (143),  (144),  and  (145)  into  Eq.  (142)  and  applying 
the  Cramdr-Rao  inequality  givea  the  reault  stated  in  the  Theorem. 
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